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Abstract 

The locomotion of Caenorhabditis elegans exhibits complex patterns. In particular, 
the worm combines mildly curved runs and sharp turns to steer its course. Both runs and 
sharp turns of various types are important components of taxis behavior. The statistics 
of sharp turns have been intensively studied. However, there have been few studies on 
runs, except for those on klinotaxis (also called weathervane mechanism), in which the 
worm gradually curves toward the direction with a high concentration of chemicals; this 
phenomenon was discovered recently. We analyzed the data of runs by excluding sharp 
turns. We show that the curving rate obeys long-tail distributions, which implies that large 
curving rates are relatively frequent. This result holds true for locomotion in environments 
both with and without a gradient of NaCl concentration; it is independent of klinotaxis. 
We propose a phenomenological computational model on the basis of a random walk with 
multiplicative noise. The assumption of multiplicative noise posits that the fluctuation of 
the force is proportional to the force exerted. The model reproduces the long-tail property 
present in the experimental data. 
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INTRODUCTION 

The soil nematode Caenorhabditis elegans ( C. elegans) is unique in that its 302 neurons and their network have 
been fully identified [1]. In addition, the entire genomic sequence of C. elegans has been determined [2]. On 
agar surfaces, the worm lies on one of its two sides and exhibits complex locomotion by crawling forward and 
backward. Complex locomotion is particularly evident when the worm responds to various types of sensory 
inputs. Examples of such functional responses include chemotaxis, thermotaxis, and isothermal tracking. The 
identification of specific neurons involved in these and other responses has been a topic of intensive research [3HI| . 

The complex locomotion of C. elegans, which underlies functional responses of the worm, accompanies 
undulations of the entire body with a period of 2-3 s. In general, trajectories of the worm are not straight 
but mildly curved. The worm often turns sharply. Sharp turning events are called under different names, such 
as reversals, sharp turns, and omega turns, depending on the specific manner of turning. A bout of these 
sharp turning events is generally called a pirouette, and pirouette is known to improve the efficiency of the 
worm's locomotion toward chemoattractants (chemotaxis) [5|. Computational modeling has been a useful tool 
for understanding the locomotion. Most existing computational models assume that sharp turning events, or 
pirouettes, are the main drive of the functional locomotion behavior of the worm [SHllj. 

Recently, another mechanism of chemotaxis was discovered [7]. It was called weathervane mechanism 
in [7], but more generally called klinotaxis |12] . In klinotaxis in C. elegans, an animal performing a run 
(segment of trajectory between adjacent sharp turns; sec MATERIALS AND METHODS for definition) senses 
the spatial gradient of chemical compounds in the direction lateral to its forward movements and gradually 
curves toward a higher concentration of the chemicals. This is distinct from the pirouette mechanism in which 
pirouettes occur frequently when the concentration of the attractant progressively decreases [8], which is a 
form of klinokinesis |12) . Klinotaxis suggests that modeling the trajectory of the worm by the combination of 
pirouettes and straight runs, which do not allow for klinotaxis, may be an oversimplification. There have been 
few computational studies of runs, except for some on klinotaxis [7] and isothermal tracking |13) . 

In this paper, we analyze run data of the locomotion of C. elegans in either the absence or presence 
of a chemoattractant in detail. Our main finding is that the curving angle per unit time exhibits a long-tail 
distribution (the precise definition of the curving angle is given in MATERIALS AND METHODS). This finding 
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is robust against tlic influence of the attractant, laser ablation of neurons, and bias in the curving inherent in 
each worm. Then, we propose a minimal computational model using a correlated random walk, which has 
been used for modeling the migration of insects and mammals |14H18] . In contrast to these previous cases, the 
increment of the curving angle cannot be modeled as a Gaussian distribution, which would arise as a summation 
of many independent and identically distributed noisy microscopic movements |19[|20j . A random walk with 
multiplicative random noise |21| explains the long-tail behavior in the locomotion of C. elegans. 

MATERIALS AND METHODS 

Experimental details 

We analyzed a subset of the experimental data used in ref. [7|. We briefly explain the experimental setup and 
the data (see ref. [7] for details). 

We used data from two sets of experiments. In one environment, each worm was placed on an agar plate 
without NaCl. This data set is the main focus of the present study. Data obtained from 46 worms in this 
environment were analyzed. In the other environment, we placed the worms on a grid format plate in which 
NaCl was spotted in a grid pattern and allowed to diffuse on chemotaxis agar as follows. Data obtained from 
53 worms in this environment were analyzed. 

Worms [C. elegans; Bristol strain N2) were cultured at 20°C using a standard technique. The agar plates 
contained 10 ml of 1 niM CaCl2; 1 niM MgS04; 5 mM potassium phosphate, pH 6.0, and 2% agar in a 9-cm- 
diameter plastic dish. For assays with NaCl (chemotaxis assays), the agar plates were prepared as follows: 1 

each of 200 mM NaCl in a low-salt buffer was spotted onto 12 points in a grid, which were spaced 2 cm 
apart; the positions of the centers of the spots were {x,y) = (10 mm, 10 mm), (10 mm, —10 mm), ( — 10 mm, 
10 mm), and so on (see [7]). The assay plates were left for 1 h after spotting the NaCl-containing buffer. 

In the tracking system, images were captured and analyzed to locate the centroid of the worm. Images 
were captured at intervals 0.4-0.6 s. Each worm was tracked for 20 min after being placed on the agar plate. 
Tracking data were discarded when a worm was considered to be immobile during an entire recording session. 
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Data analysis 

A sample trajectory of a worm in the two-dimensional plane is shown in Fig. [T]4.. The symbols represent the 
location at each image capture. We discard sharp turns (gray squares) defined in the following and analyze the 
rest of the data (dark circles). To this end, we define the displacement vector as the difference between the 
coordinate observed at a time point and that observed at the previous time point (Fig.[T}3). The angle from the 
X-axis to the displacement vector at time U measured clockwise is denoted as 6i. Following [5], we define the 
curving rate as 



This definition is slightly different from that in [7], in which the curving rate is defined by the derivative of 6 
with respect to the distance traveled. The determination of a curving rate involves three successive data points. 
The speed in addition to the curving rate is required to specify the locomotion of the worm. Although the speed 
of the worm varies with time, we neglect these variations, as in the case of previous modeling studies [711511^ 

We remove sharp turns as follows (Fig. [Ip). A sharp turn is defined as an event in which — Oi\ > 100° 
is satisfied. A sharp turn segment (shaded area in Fig. [T]C) is defined as the segment of the trajectory within 
0.1 mm from a sharp turn event. The length of the trajectory is measured by the accumulated length of the 
piecewise linear path connecting the two-dimensional coordinates. We remove all the sharp turn segments and 
analyze the rest of the data, which we call the run data. A run is a section of the trajectory between two 
successive sharp turn segments. A single worm generates a trajectory that generally includes multiple runs. 

Our data were recorded at irregular intervals. To avoid the technical complication of calculating the auto- 
correlation of the curving rate in this situation, we proceed as follows. We first generate a continuous trajectory 
on the two-dimensional plate from a sequence of the coordinates of the worm by using the third spline inter- 
polation (numerically implemented using GNU Scientific Library). Then, we resample discrete time series of 
the coordinates from the interpolated continuous trajectory every 0.3 s. Next, we detect sharp turns of the 
resampled trajectory and exclude the sharp turn segments. 4>s,i,j, ^s,i, and Mg denote the jth curving rate in 
the ith run of the sth worm (1 < s < 46), length of the ith run of the sth worm, and number of runs of the sth 
worm, respectively. We combine all the runs from all the worms to calculate a single autocorrelation function 
because reliably calculating the autocorrelation necessitates many data points. 
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The number of data points used for calculating the autocorrelation function with lag r is given by 

46 Ms 

s=l 1=1 

The average curving rate for the data points that are included in the calculation of the two-point correlation 
function with lag r is equal to 

46 M, max(f,..i-r,0) 
s—l i—1 j—1 

The autocorrelation function with lag r (i.e., 0.3 x r s), denoted by Rt, is defined by 



^ s=l 1=1 3 = 1 



46 Ms max(£a,i--r,0) 

^ E E E ('^-.^- - ) ('^-^^■+- - (^)r) 



46 Ms max(fs,i-r,0) 

rEE E (^.^.-(^u' 



46 Ms max(fs,i-T,0) 



(4) 



-EE E (<^-..+r - 



\ ^ S=l i=l i = l ^ ^ i5=l 4=1 j = l 

Note that Rq — X. Wc verified that the interpolation does not significantly affect the trajectory or the distribu- 
tions of (/). We employed the spline interpolation only for calculating the autocorrelation function. 

When NaCl is present, the concentration of NaCl created by a single spot of NaCl was evaluated by the 
solution of the diffusion equation: A^o exp(— r^/(4£'i))/(47r(i£'i), where iVg = 2.0 x 10^* (mM) is the amount 
of NaCl spotted, 13 = 1.5 x 10^'^ (mm^/s) is the diffusion constant of NaCl, d = 1.8 (mm) is the thickness of 
the plate, i is the time (in seconds) after spotting, and r is the distance (in millimeters) from the spot. The 
summation of the estimated NaCl concentration based on 12 NaCl spots defines the total concentration of NaCl, 
which we denote by C. 

To examine klinotaxis, we compute the gradient of the concentration of NaCl along two directions. One is 
the forward direction, which is parallel to the displacement vector. The corresponding gradient is denoted by 
dC/dxf. The other is the lateral direction, which is perpendicular to the forward direction. The corresponding 
gradient is denoted by ACjAxg. Wc define dC/da;£ to be positive if C on the right side of the worm, when 
viewed from above, is larger than that on the left side. 
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RESULTS 

Long-tail behavior of curving rate 

We analyze runs from two sets of experiments from different environments: ones with and without a gradient 
of the NaCl concentration (NaCl gradient for short). We mainly focus on the data recorded in the absence of 
an NaCl gradient. 

A time series of the curving rate (/) of a worm during a run is shown in Fig. cf) has a periodic component 
and accompanies large fluctuations. Periodic behavior with a period of approximately 3 s is derived from 
the body undulation. The autocorrelation function of the curving rate shown in Fig. [33 indicates oscillatory 
damping. 

The distribution of </> for a typical worm is indicated by the dashed line in Fig. The distribution has 
two peaks owing to the oscillatory undulation. The significance of the two peaks varies between worms and 
runs. When we combine the data of all the worms, the distribution of cj) is unimodal (solid line in Fig.W), 
which reflects the heterogeneity of the worms. To inspect the tails of the distribution, we divide the samples 
into a group with > and that with < 0. For a single worm, the cumulative distributions of \(f>\, defined 
in the present study as the fraction of samples among all the samples such that the absolute value is larger 
than |0|, are separately shown for the two groups in Fig. [2p. Remarkably, both two cumulative distributions 
of (f> have long tails. The long tails are also evident when we combine the data of all the worms (Fig. 03). The 
mean of for the combined data is almost equal to zero, and the standard deviation of (/) is equal to 15.8. The 
cumulative distribution of the Gaussian distribution with mean zero and standard deviation 15.8 is indicated 
by the dot-dashed line in Fig. [2p. The distribution of cf) is not accurately fitted by the Gaussian distribution 
because of the long tail present in the experimental data. A manual power-law fit to the tail of the cumulative 
distribution of \(f)\ yields the power law cx (i-e., oc |(^|~'^'^ for the original distribution). The exponent 

and the significance of the power law depend on the worm and on the sign of (jj. Nevertheless, for all the data 
examined, the tail is much longer than that of the Gaussian distribution of the same value of the standard 
deviation. 

To show that the observed long-tail distribution of (j) is not a byproduct, we perform a few tests. First, 
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the long tail is not caused by sharp turns. If we include the sharp turns in the analysis, we find a peak at a 
large value of cfi (Fig. (33) that is distinct from the long tail observed in the run data. For visibility, the original 
distributions of \(j)\, not the cumulative ones as in other similar figures, are plotted in Fig. 

Second, the long tail is not caused by the immobility of the worm. The mean of |0| is plotted against the 
instantaneous speed of the worm in Fig. |3J\. tends to be large when the worm runs slowly. Excluding the 
values of for which the instantaneous speed of the worm is small, i.e., less than 0.05 mm/s, preserves the 
long-tail distribution (Fig. [3J3). 

Third, the occurrence of a large value of \cf>\ is not restricted to periods near sharp turns. We recover 
the long-tail distribution even if we exclude a significant length of run (i.e., 6 s) before and after sharp turns 
(Fig.EP). 

In the presence of an NaCl gradient, the run is chiefly characterized by klinotaxis [7], according to which 
the worm senses the lateral gradient of NaCl and gradually curves toward the direction in which the NaCl 
concentration increases. We quantiflcd klinotaxis by following the methods developed in [7]. In Fig. the 
mean of <f) of single worms is plotted against the values of the lateral and forward gradients of NaCl concentration. 
The error bars in the figure represent the standard deviations calculated based on the mean values of from 
different worms. Klinotaxis is evident from the positive correlation between (jj and the lateral NaCl gradient 
(i.e., dC/dxi). During the run. the worm tends to curve its body to encounter a higher concentration of the 
NaCl. Note that <f> and the forward NaCl gradient (i.e., dC/dxf) are almost uncorrelated. 

Figures |4|3 and C indicate that the standard deviation and the skewncss of (/) for single worms are insensitive 
to both the lateral and forward NaCl gradients. The dependence of the skewness on the NaCl gradient (Fig.|4p), 
which was found previously [7], does exist but is negligible in the present analysis. Therefore, we expect that 
the NaCl gradient effectively shifts the distribution of cj) by an amount proportional to the lateral NaCl gradient 
but hardly changes the shape of the distribution of 0. Sample distributions of cj) conditioned by two ranges of 
dC/dxi values are compared in Fig. |3p. They do not differ considerably except for translation, as expected. 

Figure m suggests that klinotaxis translates the distribution of by a small amount and does not affect the 
long tails. A sample trajectory and the distributions of (f> obtained from the experiments in the presence of an 
NaCl gradient are shown in Fig. [5] The results are qualitatively similar to those obtained from the experiments 
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in the absence of an NaCl gradient (Fig. [2]) . 

We conclude that the long-tail behavior of the run is robustly observed both in the absence and presence 
of an NaCl gradient. 

Computational modeling 

We developed a phcnomcnological model of the run on the basis of the correlated random walk [TllSl fTl^EIIll^ . 

We indexed the discrete time by i, i.e., ti = ti-i + At, where At is a time interval. We modeled the 
reorientation of the worm at time ti by 

4>, ^ Rcos{iq^) + li, (5) 
??j = + — (ti - , (6) 

where R, rji, and Ti are the amplitude, phase, and instantaneous period of undulation, respectively. £,i is the 
noise inherent at time ti. 

If were independent Gaussian noise, as is the case for the standard correlated random walk, the long-tail 
behavior of the curving rate is not reproduced. Therefore, we assumed multiplicative noise ^i |21] . such that 

(7) 

where £ A/'(0.0, cr); is chosen from the normal distribution with mean and standard deviation a. The 
procedure for fitting the parameter values is described in Appendix I. 

The assumption of random multiplicative noise is based on the premise that a large movement accompanies 
a proportionally large noise, which is known for the arm and hand movements of humans |24H27j . To examine 
the plausibility of this assumption, we consider an alternative model in which long-tail noise is directly applied 
in each time step. If ^i in Eq. [5]is the independent long-tail noise, the long-tail distribution of (l)i is produced. 
By allowing an autoregressive (AR) term with which to calibrate the autocorrelation of 0i, we define the 
alternative model by replacing Eq. [5]by 4>i = i?cos (rji) + w(j)i-i + ^i, where < w < 1. Note that is a long- 
tail noise; the AR model with the Gaussian ^i, as is usually assumed, does not yield the long-tail distribution 
of (f>i irrespective of the parameter values and the number of autoregressive terms on the right-hand side of 
the model. We compare the two models by measuring (^{4>i — ti'(/)i_i)^^, where (•) is the average over i for a 
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2 



given |(/>i_i|. For the alternative long-tail noise model, we obtain (^{(pi ~ w(j)i-i) y = y ^i?cos (77^) + £,ij j- £,i 
is identically distributed, and its typical magnitude is much larger than that of the sinusoidal undulation (i.e., 



and serve as an estimate of the magnitude of the noise for a proper value of w. In contrast, for our multiplicative 
noise model, Eqs. [5] and [7] imply that l{4>i — increases with 



It turns out for any w that {^{(jii — w(j>i^i) y and arc positively correlated for both the experimental 

data (Fig. and the multiplicative noise model (Fig. |33). This is not the case for the long-tail noise model 
with w = 0.95 (Fig. [6j3) or any other values of w; the plots for other values of w almost completely overlap with 
those for w ~ 0.95 shown in Fig. |nj3. Therefore, at least some components of the locomotion are likely to be 
explained by the multiplicative noise. We have verified that the addition of the AR terms to the multiplicative 
noise model little affects the behavior of the model except for some changes in the autocorrelation of (pi . 

A sample trajectory obtained from the numerical simulations of the multiplicative noise model is shown in 
Fig. [3V, which is qualitatively similar to that shown in Fig.[2K. In Fig.[7j3, the autocorrelation function obtained 
from the model (dashed line labeled "constant T^") is compared with that of the experimental data (solid line). 
Even though the model reproduces the oscillatory nature of the autocorrelation function, the autocorrelation 
does not decay, indicating the excessive periodicity inherent in the model. The strong periodicity results from 
the fact that the undulation period Ti is fixed. The distribution of (j) and its long tail obtained from the 
numerical simulations roughly resemble those of the experimental data (Fig. [7p and D). 

The undulation period of the experimental data fiuctuates. To incorporate this factor, we extended the 
model by allowing Ti to fluctuate with time. We set T,; = 2.8 -I- C, where C S 7V(0. 0,0.9). We generated Ti 
until a positive value of Ti was obtained. This modification of the model yielded the exponential decay in 
the autocorrelation function, as found in the experimental data (dash-dotted line labeled "fiuctuating T^" in 
Fig.[2j3). The modified model does not quantitatively agree with the autocorrelation function of the experimental 
data, partly because the fiuctuation in Ti is crudely implemented. Note that the modification of the model does 
not change the distribution of p (dash-dotted lines in Figs.[7p and D). The sample trajectory is also similar to 
that of the original model and that of the experimental data. 

The long-tail behavior is not restricted to the specificity of the multiplicative noise, as defined by Eq. [T] 



Rcos{rii)) under the fitted parameter values. Therefore, ({(j>i — w4>i^i)^\ w (ff 




would depend little on 
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The behavior of the model is ahuost the same when we replaced Eq. [7] with a different noise model given by 

I, = e,(0.-i + (/>.-2)/2 (Fig.nD. 

Under NaCl gradients, klinotaxis seems to simply translate the distribution of (j) (Fig. 2]). Therefore, we 
extended the model such that klinotaxis adds linearly to the random walk. In other words, the curving rate 
observed under NaCl gradients is assumed to be the summation of 4'i given by Eq. [S] and the klinotaxis term 
a X dC/dxi. In the numerical simulations, we set a = 19.24, which is the value of the slope determined by the 
linear regression for the mean of the data shown in Fig. HJ^. 

The numerical results in the presence of an NaCl gradient almost overlap those from the computational 
model in the absence of an NaCl gradient shown in Fig. [T] This implies that the model reproduces the long-tail 
behavior of cf), as observed in the experimental data in the presence of an NaCl gradient (Fig. [5|). We did 
not calculate the autocorrelation function because the curving rate is nonstationary when the worm performs 
chemotaxis. We have verified that our model worm performs chemotaxis as in previous models and experimental 
data [ZIIH]. 

DISCUSSION 

We have analyzed the runs of C. elegans and discovered long-tail distributions of the curving rate. We have 
proposed a phcnomcnological model that qualitatively reproduces the behavior of the curving rate. Our model 
is a type of correlated random walk with multiplicative noise. 

The long-tail behavior is robust against various types of perturbations. First, it is observed both in the 
absence and the presence of NaCl gradients such that it is independent of klinotaxis. Second, our results are not 
restricted to the specific version of the multiplicative noise model assumed in this study (see Fig.|8]). Third, the 
long-tail behavior is robust against laser ablation of various neurons. As shown in Fig. |9l the long-tail behavior 
is preserved regardless of whether the ablated neurons are responsible for klinotaxis for NaCl chemotaxis (AIZ 
neurons; see [7]) or not (RIA neurons). This result gives another support to the observation that the long-tail 
behavior is independent of klinotaxis. We speculate that the long-tail behavior reflects intrinsic properties of the 
neural circuit for head oscillation, which is currently unidentified. Fourth, the bias in the curving rate that each 
worm owns jTKH], perhaps because the worm crawls on one of its sides, does not affect the long-tail behavior (see 
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Appendix II) . We have also verified by numerical simulations that the long-tail behavior of 4> does not hamper 
klinotaxis. 

The travel length per unit time obeys the Levy distribution (i.e., a power-law distribution) for various 
animals. Such a Levy flight is suggested to enhance the efhciency of foraging PSII^ . The long-tail distributions 
in the curving rate that we have discovered arc not directly related to the Levy flight. In the Levy flight, the 
distance of the movement per unit time follows a long-tail distribution. Accordingly, long jumps occasionally 
occur to possibly accelerate the foraging. In the locomotion of C. elegans, the distance of the movement per unit 
time is narrowly distributed, which reflects the fact that the worm moves only by means of regular undulation 
of the body. A long-tail distribution of the curving rate docs not appear to affect the long-term locomotion 
behavior relevant in foraging or taxis. 

The trajectories of the worms are highly noisy. The correlated random walk is a useful tool for modeling 
noisy migrations of C. elegans [711511^ and other animals P^lSlllSfBO] . The noise per unit time is often modeled 
as Gaussian noise [Zllin]- However, Gaussian noise does not lead to distributions with long tails. Obviously, 
we could reproduce the long tail if independent noise obeying the long-tail distribution is applied at each time 
step |17j . However, this independent- noise model contradicts the positive correlation between the statistics of 
the consecutive curving rates present in the data (Fig. [6]). To be consistent with both the long-tail behavior and 
the positive correlation, we have assumed that the amplitude of the noise is proportional to the curving rate. 
This assumption is based on the premise that exerting a large force accompanies a proportionally large noise. 
This is the case in the motor activities of human beings; the amplitude of the physiological noise is proportional 
to the magnitude of the hand and arm movements or the generated force [24H27] . We note that the positive 
correlation in the data is not as strong as the expectation from a random walk with purely multiplicative noise. 
The locomotion of the actual worm may be generated by combination of multiplicative noise and independent 
noise. To the best of our knowledge, the relationship between the magnitude of the force and that of the noise 
has not been analyzed for animals other than humans. An investigation of this issue may be an interesting 
future problem. Explicit modeling of the animal's movement should also be attempted in the future. 
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Appendix I: Determination of parameter values of the 
computational model 

We determined the four fitting parameters of the model, i.e., At, R, Ti, and a, as follows. On the basis of 
the long-term oscillatory period of the data, we set Ti = 2.8 s independent of i (also see Fig. UK). We also 
investigated the case of fluctuating Ti in the main text. The amplitude of the oscillatory head movement is set 
to i? = 10.0°/s so that the dual peaks in Fig. [5p are roughly recovered. These two parameters are necessary 
for reproducing the undulation of the worm. The standard deviation of the intrinsic noise (i.e., a) controls 
the dispersion of (f>. In order to match the standard deviation of the distribution of (j) with the experimental 
one, we set cr = 0.9. The parameters R and a mainly control the distribution of (j) for small and large \(f)\, 
respectively. It was found that a does not affect the timescale of the autocorrelation to a great extent. At is 
a relevant parameter of the present model and mainly affects the timescale of the decay of the autocorrelation 
function of (j) for small time lags (see Fig. [TU)) . To fit the experimental data in this aspect, we set At — 0.6 s. 

Appendix II: Worm-dependent bias in curving rate 

The standard deviation of (f> in the absence of an NaCl gradient is plotted against the mean of cj) in Fig.[TTJ One 
data point corresponds to one worm. Figure [TT1 indicates that the bias (i.e., mean (jj) inherent in each worm is 
small relative to its standard deviation. However, the effect of this bias is sufficiently large to generate a visible 
rotating trend of the locomotion. 

To examine the effect of the bias on our main results, we considered the curving rate generated by adding 
a constant bias as well as the klinotaxis term a x dC / dxg to the right-hand side of Eq. [5] [HIUS] . We confirmed 
that this modification of the model does not significantly affect the distribution of and the shape of the 
autocorrelation function, regardless of the magnitude of the constant bias. The bias has a significant effect on 
a long time scale, but not on a short one. Because our analysis of (j) is concerned with a short time scale, we 
did not detrend the bias in the analysis of our experimental data in the main text. 

Note that the concept of the bias discussed here is different from the tendency that the worm migrates in 
a specific direction in which the concentration of the attractant increases |18] . 
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FIGURE LEGENDS 

Figure 1: (A) Sample trajectory of the worm. (B) Definition of the curving rate (pi. (C) Definition of the sharp 
turn segment (shaded area). 

Figure 2: Resuhs obtained from experiments in the absence of an NaCl gradient. (A) Sample trajectory of 0. 
(B) Autocorrelation function of (p. (C) Distributions of cj) based on the runs of a single worm and those of all 
the worms. We use a bin width of 4°/s . (D) Cumulative distributions of \4>\ (i.e., fraction of values above 101) 
on the log-log scale. The results for one worm and the combined one for all the worms are shown separately 
for positive and negative values of (f>. The dot-dashed line in (D) indicates the cumulative distribution of the 
Gaussian distribution with mean zero and standard deviation 15.8. (E) Distribution of |0| for the data before 
and after sharp turns are excluded. The data for all the worms are combined to generate each distribution. 

Figure 3: (A) Relationship between the mean of \4)\ and the instantaneous speed of the worms. The mean 
values are calculated from all the runs of a single worm by categorizing the data into bins of width 0.025 mm/s 
and taking the average in each bin. Then, the mean values from different worms for the same bin are used for 
calculating the grand mean and standard deviation. The grand mean and error bar defined by the one standard 
deviation range are plotted in (A). (B) Cumulative distribution of |0| when data points for which the speed 
of the worms is less than 0.05 mm/s are excluded. (C) Cumulative distributions of |0| when the data points 
within 6 s before or after each sharp turn are excluded. For each figure, the data obtained from all the worms 
are combined to generate the cumulative distribution. 

Figure 4: Klinotaxis. (A) Mean, (B) standard deviation, and (C) skewness of (p for different values of lateral and 
forward NaCl gradients. The mean, standard deviation, and skewness for a value of the concentration gradient 
are calculated from all the runs of a single worm. Then, the values from different worms for the same value of the 
concentration gradient are used for calculating their mean and standard deviation. The mean values, together 
with the error bar defined by the one standard deviation range, are plotted in (A)~(C). (D) Distributions of (p 
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conditioned by two ranges of lateral NaCl gradient dC/dxf,. We combined the data of all the worms to make 
the two distributions sufficiently smooth. 

Figure 5: Results obtained from experiments in the presence of an NaCl gradient. (A) Sample trajectory of (j>. 
(B) Distributions of (p. (C) Cumulative distributions of \(f>\. See the caption of Fig. [2]for legends. 



Figure 6: Testing the multiplicative noise hypothesis. We plot y ~ w(f)i-i) j against l^ii-il for four values 
of w for (A) experimental data and (B) random walk model with multiplicative noise. The results for the AR 
model with w = 0.95 and the power-law noise distribution with exponent 3 arc also shown in (B). Each plotted 
point represents the square root of the average of {cj)i — wcj)i-if' in a bin. Each bin has width 20°/s and contains 
the corresponding values of |<?I>i-i| obtained from all the worms. 

Figure 7: Results obtained from numerical simulations of the computational model in the absence of an NaCl 
gradient. (A) Sample trajectory of when Ti is constant. (B) Autocorrelation function of (p. (C) Distributions 
of (j). (D) Cumulative distributions of |(/)| on the log- log scale. In (B), (C), and (D), the experimental data, the 
numerical results with constant T^, and those with fluctuating Ti arc compared. 

Figure 8: Numerical results in the absence of an NaCl gradient when a different noise model is used. Eq. [7]is 
replaced by = ^i((/>i-i -t-0i_2)/2 with G 7V(0.0, 1.24). We modified the standard deviation of the noise from 
that of the original model (Eq. [7|) to obtain a reasonable fit to the experimental data. Long-tail distributions 
of 4) arc produced, as shown in D. See the caption of Fig. [7] for legends. 

Figure 9: Cumulative distribution of |0| for the worms subjected to laser ablation of either of two different 
neurons (RIA and AIZ). The laser-ablated worms are subjected to an NaCl gradient. Each distribution is 
generated by combining the data of all the worms. The long-tail behavior of the run is supported by these data. 
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The scaling exponent of the long-tail distribution depends on the type of ablation. The results in the absence 
of ablation are also plotted for comparison. 

Figure 10: Numerical results for different values of At. (A) Cumulative distributions of \<f)\. The distribution 
of the curving rate is nearly independent of the value of At in the computational model. (B) Autocorrelation 
function of (j). The autocorrelation function for small time lags depends on At. A small At yields a rapid decay 
of the autocorrelation function for a small lag. In the main text, we set At = 0.6 s such that the autocorrelation 
function for small time lags is close to that of the experimental data. 

Figure 11: Relationship between the standard deviation of (f) and the mean of (f) for each worm in the absence 
of an NaCl gradient. 
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